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CURRENT RESERVOIRS IN THE SIMPLE EXCLUSION PROCESS 

A. DE MASI, E. PRESUTTI, D. TSAGKAROGIANNIS, AND M.E. VARES 

Abstract. We consider the symmetric simple exclusion process in the interval [—A, A] 
with additional birth and death processes respectively on (N—K, A], K > 0, and [—A, — A+ 
• K). The exclusion is speeded up by a factor A 2 , births and deaths by a factor A. Assuming 

propagation of chaos (a property proved in a companion paper, [3]) we prove convergence 
in the limit A — > oo to the linear heat equation with Dirichlct condition on the boundaries; 
' the boundary conditions however are not known a priori, they are obtained by solving a 

non linear equation. The model simulates mass transport with current reservoirs at the 
boundaries and the Fourier law is proved to hold. 

a,: 

d ' Contents 



> 



^1- 



1. Introduction 1 

2. Model and main results 2 

3. The discretized evolution 6 



1 4. Probability estimates for a random walk with reflections 7 



5. Proof of the hydrodynamic limit 10 

C<S : 6. Proof of the Fourier law 13 

' References 19 

O ' 



1. Introduction 



Stationary non equilibrium states are characterized by the presence of steady currents flowing 
through the system and a basic question in statistical mechanics is to understand their 
structure. Many papers have been devoted to the subject, we just quote a few of them 
where the issue is addressed in the context of stochastic interacting particle systems, [1], 
[2], [4], [5]. We want to produce currents by acting only at the boundaries [of the region 
where our particle system is confined], so that the bulk dynamics is left unchanged. There 
are two natural ways to proceed. Thinking of a one dimensional system in an interval 
(this paper is indeed about the one dimensional symmetric simple exclusion process, SSEP) 
we may produce a current by sending in particles from the right at some small rate (i.e. 
inversely proportional to the size of the region) and taking out at same rate particles from 
the left. Alternatively we may add and subtract at unit rate (i.e. much faster than previously) 
particles on the right trying to keep fixed a given density p + close to the right boundary; 
same is done on the left and if the two densities p± are different, say p + > p_, then we have 

l 
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a positive density gradient which by the Fourier law induces a negative current (inversely 
proportional to the size of the region). 

While the second mechanism has been much studied in the context of the Fourier law and to 
investigate the invariant measures when stationary currents are present, as in the references 
quoted above, to our knowledge the first mechanism has not been examined so far, even 
though it may look the most natural and direct to produce a current and our purpose here 
is to start its analysis in the simplest possible context. For this reason we consider the 
d — 1 SSEP in an interval Ajv = [— N, N], N a positive integer (we are interested in the 
asymptotics as N — > oo). The process takes place in {0, 1} Ajv (at most one particle per 
site) and time is speeded up by a factor N 2 (to match the length of the interval Ajv). Thus, 
independently each particle tries to jump at rate N 2 /2 to each one of its n.n. sites, the jump 
then takes place if and only if the chosen site is empty, see the next section for a formal 
definition; jumps outside An are suppressed. To induce a current we modify the process by 
sending in from the right and taking out from the left particles at rate Nj/2, j > a fixed 
parameter independent of N (to compare with the statements in the beginning we should 
divide the rate by A^ 2 because times have been speeded up by a factor A^ 2 ). As we want the 
boundary processes localized at the boundaries we fix two intervals I± of length K at the 
boundaries and we send in particles only in I + and take out particles only from 7_. It may 
however happen that I + is already full or J_ empty, then our mechanisms abort, so that the 
current really flowing in the system will not be exactly what desired (but hopefully close if 
K is large). This seems unavoidable if we insist to localize at the boundaries the birth-death 
processes or to consider lattice gases. In this paper K may be arbitrarily large but fixed; we 
cannot allow K to grow with N and a system of positive and negative charges (for instance 
a zero range process with positive and negative particles) requires a different approach. 
In this paper we derive the hydrodynamic equations for our model (in the limit as N — > oo) 
under the hypothesis of "propagation of chaos" , a property proved in a companion paper, [3]. 
The hydrodynamic equation is the linear heat equation in the "macroscopic" interval (—1, 1) 
with Dirichlet boundary conditions at time and at ±1: the values at ±1 are not a-priori 
given, they are obtained by solving a coupled system of two non linear integral equations. 
We also prove the validity of the Fourier law; in particular the currents which enter and exit 
from the system are at all times equal to the local density gradient at ±1. 
We hope to continue this research project by studying the fluctuations field in the hydrody- 
namic limit and then the structure of the invariant measure in the limit as A" — > oo. 

2. Model and main results 

Notation and definitions. 

A A? is the interval in Z with endpoints ±AT, denoted by A/y := [-N, N]. We write e = 1/N, 
fix an integer K > 1, write 1+ = [N - K + 1, N] and J_ = [-N, -N + K - 1]. Particle 
configurations are elements r\ of {0, 1} Ajv , r\{x) = 0, 1 being the occupation number at x G A/y. 
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We shall study the Markov process on {0, 1} Ajv with generator L e := e~ 2 (L + eL b ), where 
L b = L b:+ + Lb - and 

N-l 

Lof(v) ■= j E lf(v (x ' x+1) ) ~ f(v)\, Wfo) := 3 -J2 D ±V (x)[f( V M) - f( V )\ , (2.1) 

x=-N x£l± 

T)( x > being the configuration obtained from 77 by changing the occupation number at x, rj^ x,x+1 ' 
by exchanging the occupation numbers at x, x + 1; for any u : A N — > [0, 1] 

D + u(x) = [1 - u(x)]u(x + l)u(x + 2) . . . u(N), x e 1+ 

D^u(x) = u(x)[l-u(x-l)][l-u(x-2)]...[l-u(-N)], 16L (2.2) 

L is the generator of the SSEP (and of the stirring process as well). L 6)+ and L bj _ are 
generators of birth respectively death processes, the former is active in I + the latter in J_. 
The parabolic nature of the stirring process suggests to scale time as the square of space, 
hence the factor e~ 2 in the definition of L e . It readily follows from the structure of the 
generators that the expectations K e [r](x, t)] satisfy the relations 

^E e [r](x,t)} = ^A e EM^t)] + e- lJ -(l xeI+ E e [D +V (x,t)] - l xeI _E e [D. V (x, t)] ) (2.3) 

where A e = e~ 2 A, A the discrete Laplacian in with reflecting boundary conditions: 

Au(x) = u(x + 1) + u(x — 1) — 2m(x), \x\ < e _1 
Au(±N) = u(±(N - 1), t) - u(±N, t). (2.4) 



Propagation of chaos. 

Due to the last term, (2.3) is not a closed equation in E, e [r)(x, t)], but since the stirring 
generator is the leading term in L e and the invariant measures for the stirring process on the 
line are product Bernoulli measures, it looks natural to conjecture "propagation of chaos" , 
i.e. that the measures at time t > are approximately product (as e — > 0). If the law at 
time t > were a true product measure, then, instead of (2.3) the expectations ~E e [r](x, t)] 
would satisfy the closed equation: 

^Pe(x, t) = ^ A e p e (x, t) + e~ l3 - (l xeI+ D + p e (x, t) - l x£l _D_p € (x, t)j , (2.5) 
which will be referred to as "the discretized time evolution" . 

The Cauchy problem for (2.5) with [0, 1] -valued initial datum p e (-,0) has a unique global 
solution which also takes values in [0, 1]. Indeed (2.5) is a first order system of ordinary 
differential equations with polynomial non linearity hence local existence and uniqueness. 
Global existence follows because the solution has values in [0, 1], which in turns is a conse- 
quence of the fact that D±u(x) vanishes when u(x) = 1, respectively u(x) = 0. A formal 
proof is given in Proposition 3.1. 
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Hydrodynamic limit. 

The first result in this paper (proved in Section 5) shows that p e converges as e — > to a 
limit function which then identifies the hydrodynamics of the system. 

Theorem 2.1. Suppose that p € (x,0), x G A^r, with values in [0,1], converges weakly as 
e to p(r, 0) G L°°([-l, 1], [0, 1]) in the sense that 

lime V p e (x,O)0(ex) = / p(r, 0)<f>(r)dr, for any G L°°([-l, 1], R). (2.6) 

T/ien t/iere zs p(r,t), r G [—1, 1], £ > so that for any t± > to > 0: 

lim sup sup \p e (x, t) — p(ex,t)\ = 0. (2.7) 
xeA N t <t<t! 

The function p(r, t) solves and is the unique solution of the integral equation 

p(r,t)= [ P t (r } r')p(r' : 0)dr' + J - [ |p s (r, 1)(1 - p(l, t - s) K ) 
'< 1,1] 1 Jo k 

-P s (r, -1)(1 - (1 - p{-l,t- s)) K )}ds, (2.8) 



where P t (r,r') is the density kernel of the semigroup (also denoted as P t ) with generator A/2, 
A the laplacian in [—1, 1] with reflecting, Neumann, boundary conditions, (see the Remarks 
below). 



Remarks. 

• The density kernel P t (r,r') can be expressed in terms of the Gaussian kernel 

e -(r-r') 2 /(2t) 

G t (r, r') = = — , r, r' G R, (2.9) 

V 27rt 

as follows: if if} : R — > [—1, 1] denotes the usual reflection map, i.e. ip(x) = x for x G [—1, 1], 
ip{x) = 2 — x for x G [1,3], with if; extended to the whole line as periodic of period 4, then 

P t (r,r') = G t(r,r") for r' ^ ±1 (2.10) 

r":ip(r")=r> 

P f (r,±l) = 2G t(r,r"). 

r":ip(r")=±l 

• From the expressions above and (2.8) it follows that p(-,t) is "smooth" for any t > 0: we 
are calling "smooth" a function f(r), r G [—1, 1], if it is C°° in (—1, 1), continuous in [—1, 1] 

d n f(r) 

and if for each n exist the limits — ; as r — > ±1. 

dr n 
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• Since p is smooth we can write (2.8) in differential form: it then becomes the heat equation 
with Dirichlet boundary conditions: 

§£P(r, t) = l§^p(r, t), re (-1, 1), t > (2.11) 
p{r,0) = u {r), p{±l,t) = u ± {t) 

However the boundary conditions u±(t) are not a-priori known, they must be obtained by 
solving a nonlinear system of two integral equations: 

u±(t) = [ {p(s)f±(u ± (t - s)) - q(s)f T (u T (t - s))}ds + w±, t (2.12) 
Jo 

/ + («) = f(l-«*), /_(«) = |(l-(l -«)*), 
where, writing G(r) = G(0,r), the latter as in (2.9), 
p(t) = 2^G t (4fc) ) q(t) = 2j2Gt(4k + 2) 

k&L k&Z 

(2.13) 

w +' t = Yl [ u o( r ') 2G t(l -r' + 4k)dr', w- 1 = Yl f u Q (r')2G t (r' + I + Ak)dr' 

• By a simple computation one can check that 

= j{l - p(l, t)«), = j(l - (1 - p(-l, t)f). (2.14) 

This remark will be important in the analysis of the Fourier law. 

• To characterize the asymptotic behavior of the invariant measure as N — > oo it will be 
important to study the evolution starting from arbitrary initial configurations 77^. Since 
the functions f^(r) := ?7^ Ar ^([e~ 1 r]) are in a ball of L 2 ([— 1, 1], M) they converge weakly 
by subsequences to an element of L 2 ([— 1, 1],M) and we can then apply to any convergent 
subsequence Theorem 2.1. 

• The identification of (2.8) and (2.11) as the hydrodynamic equation of the system is based 
on the assumption that p e gives an accurate description of the process. This is indeed correct 
because the "empirical averages" are close to the functions p e in the following sense. There 
is r > so that calling Jm(x) = [x — M, x + M\ fl A N , M the integer part of N a , a G (0, 1), 
then for any to > 

limlim sup supP e sup | 1 ^ My,*) - p e (y, t)}\ > 5 =0 (2.15) 



5 ~" 0e ^ t Q <t<T\oge- 1 V L x€A N \JM 



Jm(x), 



(2.15) follows from Theorem 6.1 using the Chebishev inequality, Theorem 6.1 is proved in a 
companion paper, [3]. 
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Fourier law. 

In the "hydro dynamic limit literature" the limit function p(r, t) is usually interpreted as the 
"density profile" at time t: this comes from attributing to each particle a mass e so that 
K t [erj(x,t)} is the average mass in the interval [x — |,x + |] which in macroscopic units has 
length e (as [-N, N] in the macroscopic limit shrinks to [—1, 1]). Thus K e [t](x, t)} is the mass 
density, which in the limit converges to p(r,t) (when ex — > r). Analogously, the expected 
current through a point is the average signed mass crossing that point per unit time. Let x 
be away from the boundaries in the sense that \x\ < N — K. Then it follows from (2.1) that 
the expected current through x + | is 

j (e) (*, t) = ^ E e [e{ v (x, t) - nix + 1, *)}] = -i Ee [ rts + Mj-qOM) ] (2 16) 

By a similar argument the expected currents through N and — N are: 

= X>[eD ± fo(y,f))] = -\ J2^[ D ^(y,t))]- (2-17) 

By (2.16) j^ € \x, t), \x\ < N — K, is equal to — | times the discrete gradient of the density 
in agreement with the Fourier's law, which is then satisfied before the macroscopic limit 
e — > 0, (but not necessarily in the limit, as this requires that the limit of the derivative is the 
derivative of the limit). One would expect that also j± (t) are equal to — | times the discrete 
gradient of the density, at least in the limit as e — > 0. This is settled in the next theorem 
where using the factorization properties proved in [3] we show that the limit of the current 
is both in the bulk and at the boundaries equal to — | times the gradient of the density. 

Theorem 2.2 (Validity of the Fourier law). Suppose that the process starts with a product 
measure p e such that p e [r](x) = 1] = p(ex, 0) with p(-,0) G C([— 1, 1], [0, 1]) with bounded 
derivative in (—1,1). Let p e (x,t) be the solution of (2.5) starting from p e (x,0) = p(ex,Q). 
Then, for any t > and r G (—1, 1), denoting by [u] the integer part of u, we have 

limj^([ e - 1 r],t) = -i%^. (2.18) 
2 ar 

Moreover for any t > 

limj£ e) (t) = -i^ • (2.19) 

e^o J± w 2 dr r=±l 

Theorem 2.2 is proved in Section 6. 

3. The discretized evolution 
We begin the analysis of (2.5) by proving: 
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Proposition 3.1. The Cauchy problem for (2.5) with initial datum p e (x, 0) G [0,1] has a 
unique global solution p e (x,t). Moreover p e (x,t) takes values [0, 1]. 



Proof. Write 



D*,u(x) = (1 - u(x)) u(x + l)u(x + 2) . ..u(N) 



x e 1 4 



D*_u{x) = u(x) (1 - u(x - 1))(1 - u(x - 2)) ... (1 - w(-iV)) 



x e L 



If < w(x) < 1 then D± = D±. A local existence and uniqueness theorem holds for the 
Cauchy problem (2.5) as well as for the problem with D± replaced by D± (as these are 
Lipschitz functions of u in the sup-norm topology). Denote the solution of the latter by 
p*(x,t), t < t, r > 0, recalling that the initial datum p e (x, 0) verifies < p e (x, 0) < 1 for 
any i 6 Ajy. We shall next prove that < p*(x,t) < 1 for all x and t < r. Define u(s) = 
max^gA^ p*(x, s) and suppose by contradiction that there is T < r such that u(T) > 1. Then 
there is t < T so that (i) u(t) > 1 and (ii) du{t)/dt > (because u(0) < 1 and it cannot be 
that du(s)/ds < for almost all s < T such that u(s) > 1). Moreover there exists x such 
that (a) p*(x,t) = u(t) and (b) du(t)/dt = dp*(x,t)/dt. All that leads to a contradiction 
because dp*(x,t)/dt = |A e p*(x, t) + {{D* + - D*_)p*{x, t) < 0. Indeed A e p*(x,t) < 0, because 
(x,t) maximizes p*(-,t). D* + p* e {x,t) = if x ^ I + and < in I + , because p*(x,t) > 1. 
D*„p* = if x <£ J_ and > in J_, because p*(x,t) > 0. Thus (D* + - D*_)p* e (x,t) < 0. 
Analogous arguments show that the solution cannot exit [0,1] through 0, hence p*(x,t) G 
[0,1]. As a consequence D*±p* e = D±p* and therefore p* solves (2.5) as well. By iteration, 
the previous argument extends to all times. □ 

We shall study (2.5) in its integral form: 

ds(Y,P^&y) D +Pe(y,t-s) 

-Y,Pl e Xx,y)D-Pe{y,t-s)) (3.1) 

yel- 

where is the semigroup with generator |A e , P^ e \x,y) its kernel: 

pW :=e iA. tj P^\ X ,y)=P^\y, X ) (3.2) 

The analysis of (3.1) will exploit the nice regularity properties of P t (x, y) which are estab- 
lished in the next section. 



4. Probability estimates for a random walk with reflections 

In this section we shall consider a simple random walk on A^v which jumps with intensity 
e~ 2 /2 to each of its n.n. sites, the jumps outside A N being suppressed. We denote by P^ 
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its law and call Ql the law of the corresponding unrestricted random walk on the whole Z. 
In the sequel we shall prove (in many cases just recall) bounds and estimates on P^ which 
will be used in the next sections to prove Theorems 2.1 and 2.2. We start by relating 
and Q{ , through a "reflection map" from Z to A^r which is a discrete analogue of the map 
■0 defined in the first remark after Theorem 2.1. Since the jump rate from ±N is e _2 /2 for 

P^ and e~ 2 for to relate the two it just suffices to identify N + 1 with N (as well as 
— N — 1 with —N), then the jumps of from N to N + 1 and —N to —N — 1 are like 
suppressed. We thus define: 

Definition. The "reflection map" ip^ '■ Z — > A^r is: 

• |x| < N: ipN^x) = x 

• x < —N: ipN^x) = —tPn(—x) 

• x > N: ip N (N + k) = N - (k - 1), for k = 1, . . . , 2N + 1, 
ip N ((N + 2N+ 1) + k) = -N+ (k - 1)), for k = 1, . . . ,2N + 1, 

ij N ((N + 2(2N + 1) + A;) — N — (k - 1)), for jfe = 1, . . . , 2N + 1 and so on. 

Proposition 4.1. With the above notation, 

y-ipN(y)=z 

Proof. Let / be ipN measurable, i.e. /(^at(^)) = f(x), and g the function on A^ defined 
by setting g(x) = f{x), x G Ajv. Calling Lq and Lp the generators of Qt and Pf we have 

L Qf{%) = L P g(tp N (x)) 

so that e LQt f = e Lpt g, hence (4.1). □ 
By the local central limit theorem (see for instance [6]): 

Theorem 4.2. There exist positive constants Ci, ...,c^ so that 

\Qf{x,y)-G € - H {x,y)\ < -p=G e - H {x,y), \x-y\< (e" 2 t) 5 / 8 , 

(4.2) 

Q { t\x,y) < min|c 2 e- C3|a: - J/|2/(e " 2t) ,C4e- |?/ - a;|(log|j/ - a;| - C5) }, \x - y\ > (e~ 2 t) 5/8 , 
G t being the Gaussian kernel defined in (2.9). 



The next corollary follows directly from Theorem 4.2 and Proposition 4.1. 
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Corollary 4.3. For any T > there exist c so that the following holds. 

• For all e, all t £ (0, T] and all x, y in An, 

P} e \x,y)<cG e -2 t (x,y). 

• For all e, all t £ (0, T] and all -N < x < N — I, 



G e -n(x,y). 



P t {e \x,y)-P t ie) (x + l,y) 
For all e, all t £ (0, T] , s > and all x £ Ajv, 

E |^tS(^2/) - ^ (e) (^,2/)| < 



>(0 



P t ^(x,y)-P, (e) (x,y) 



< c 



V 7 



" 2 t ' 



(4.3) 
(4.4) 

(4.5) 
(4.6) 



Proof. (4.3) and (4.4) follow directly from (4.2). By (4.3) and (4.4) we can bound the left 
hand side of (4.5) by 



E E P . C 'W) Pi 6 \z,y)-P t ie \x,y) < c E G 



\z — x\ 



-2 S X, Z 



hence (4.5). (4.6) is obtained similarly, recalling that G t (x,y) < ct~ l l 2 . 



□ 



In the proof of Theorem 2.1 we shall use the following convergence results: 



Lemma 4.4. As in Theorem 2.1 suppose that p e (-,0) converges weakly to p(-,0). Then, for 
any t > 



lim sup I £P<W(,,0)-/ P, 



^ex, r)p(r, 0)<ir 



0. 



(4.7) 



Proof. By (4.4) the family of functions f e (r): 

/ e (r) := E P} 6 \[e- 1 r],y) P e(y,0), r £ [-1,1] 



J/GAjv 



is uniformly Lipschitz, it will therefore suffice to prove pointwise convergence. We thus fix 
r* £ [-1, 1] and take x = [e~V]. By (4.1) and (4.2) 

E Pt {e \x,y)Pe(y,0) = E P*(V>°) E Qt\*>*) 

y&A N yeAjv z-.ip N {z)=y 



Epe(y,0) E G e - H (x,z) + 

J/SAjv z:ip N (z)=y 



0([e-Hy l ' 2 )+0(e- 
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Call the discrete analogue of the reflection map ip of the Remarks after Theorem 2.1, 
i.e. ^>n(x) = Ntp(x/N) for x G Ajy. It differs from tp^ by shifts and we have: 



^2 G e -2 t (x,z)- ^ G e -2 t (x,z) 

z:ij> N (z)=y z:^ N (z)=y 



< 



But (see (2.10)) 



^2 G e -2 t (x,z)= ^ eG t (ex,r') = eP t (ex,ey) 

z.ty jv {z)=y r':tp(r')=ey 



and we conclude that 



< 



| J2 Pt ie \x,y)p e (y,0)-eJ2 Pt(ex, ey)p e (y,0) 
yeAjv y&^N 

Let x = [e _1 r*], then by (2.6) with 0(r) := P t (r*,r) we have 

lim V P t {t \[e- 1 r*],y)p e (y,0) = f P t (r* ,r)p(r,0)dr (4.9) 

which proves pointwise convergence and hence the lemma as argued at the beginning of the 
proof. □ 

Lemma 4.5. Let h e (t) be a continuous function with values in [0, 1] which converges point- 
wise to h(t). Then for any t > 0, r G [—1, 1] and y G i+ 



lim / e _1 P s (e) ([e -1 r], ±y)h e (t - s)ds = f P s {r, ±lWt - s)ds. 
e ^°Jo Jo 



(4.10) 



Proof. Again, this follows easily from (4.1) and (4.2), after recalling also (2.10). Details are 
omitted. □ 



5. Proof of the hydrodynamic limit 

In this section we shall prove Theorem 2.1. We start by proving equicontinuity which is a 
direct consequence of the estimates of the previous section: 

Proposition 5.1. For any T > there is a constant c so that for any solution p e (x,t) of 
(2.5) with p e (-,0) G [0, 1] the following holds. For any x G [—N,N — 1], any t G (0, T\ and 
any e > 

\p e (x,t) -p e (x + l,t)\ <min|l,c(elog + (e- 2 t) + - ? i=)} (5.1) 
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where log + u = maxjlogw, 1}. For any < s < t, x G and e > 0: 



\p e (x, t) - p e (x, t + s)\ < min |l, c(^? + Vs lo s(^)) } 

Proof. By (4.4) 

\J2(Pt (£ \x,y)-Pt (e \x + hy))Pe(y,o)\< 



(5.2) 



y 

and for any y G i+ U 7_, 



/ e - 1 |pW(x,y)-Pi £ )(x + l,y)|^<e+ / 



hence (5.1). By (4.5) 



E ( P &(x,y) - P t i£) (x,y))p e (y,0)\ < 



By (4.6) for any y G i+ U J_ and denoting by /(£) := D±p t (y,t), y G 7±, 



e" 1 



S P s ( , e) (x, y)/(t + s - s')ds' - f Pj, e) (x, y)f(t - s')ds> 



S' \ ' 9/J V 1 ~ . s 

^0 



< e' 1 ( jf * P«(x, y)^ + J* P^(x, y)ds> + jf |P s ( , e | s (x, y) - pM(x, y)|d/ 

<c(v^+v^log(^)) (5.3) 

□ 

We now turn to the proof of Theorem 2.1, we fix T > and study the evolution in the finite 
time interval [0,T]. Since we only have that p e (-,0) —> p(-,0) weakly, it is convenient to 
introduce a regularized equation. We denote by p e (x, t\u, s), t > s > 0, the solution of (2.5) 
for t > s with u the initial datum at time s, u = u(x), x G A N ,u(x) G [0, 1]. With such 
notation we then set for any 5 G (0, T) 

Pe (M) -\*MA-.*).*) (5 - 4) 

By Proposition 5.1 the family of functions (r, £) — > ([e _1 r], t), r G [— 1,1], t G [5, T] is 
equicontinuous and bounded, hence it converges in sup norm by subsequences to a limit 
function it"(r,t). By Lemma 4.4 

u {s) {r,5)=[ Ps(r,r')p(r',0)dr'. (5.5) 
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Moreover for any integer < m < K 



lim sup 



8<t<T 



D+pW (N -m,t)-(l- u^(l, t))u^(l,t) m } 







(5.6) 



lim sup 



8<t<T 



D^pi 5 \-N + m,t)-u {5 \-l,t)(l-u^(-l,t)) m } 



0. 



By Lemma 4.4 and Lemma 4.5 it then follows that 

u (*>(r,i) = u(r,t\u^\-,8),8), (5.7) 

where the latter is the solution of the limit equation in the time interval [6, T] with initial 
datum at time 5 equal to u^(-,5). Uniqueness can be easily proved, it follows also from 
(5.9) below where we prove that the solution depends continuously on the initial datum. By 
uniqueness it then follows that p{ converges in sup norm to as e — > and not only by 
subsequences. 

We shall next study the dependence on 5 and define for t G [S, T] 

hW(t) = sup \p^(x,t) - p e (x,t)\, h {5 \t) := sup |u ((5) (r,t) -u(r,t)\ (5.8) 

zGAjv r|<l 

We are going to prove that there is ct so that for all e and 6 positive 

hi s \t) + h {s \t) <c T V6 (5.9) 

(which in particular implies uniqueness of the solution of (2.8)). It follows from (5.9) that 

limsup sup sup \p e (x,t) — u(ex, t)\ < c\/S (5.10) 
e->o te[s,T] xeA N 

which then proves Theorem 2.1. 

Proof of (5.9). From (4.3) it follows that he S \5) < cy/6, then using again (4.3), 

rt i 

h^(t)<cV5 + C -^h { P{t- s)ds. (5.11) 
Js V s 



By iteration, 



oo 

hi 5) (t) < c\f8~(l + J2C n a n {t - 5)) 



n=l 

t j pt-Si pt—Sl S n -l J 



a n (t) := / —=ds 1 / —=ds 2 ... I —=ds n 
Iq v s i Jo v s 2 Jo V s " 



By (5.12) below we then have hi 5 \t) < d\f8\l + e 77 *" 21 ^. Same argument applies to h^(t) 
and (5.9) is proved. 
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Lemma 5.2. With a n (t) as above, 



a n (t) < (Trt)Se-t [lo ^- 1] 



(5.12) 



Proof. We have 

f n 1 

a n{t) = I l S iH hsn<£ I I <^ s l • • • ^ s n- 

•/[0,t]» ~ = i V s * 

We change variables by setting Si — tit and get 

f - 1 

a„(t) = (Vf) n / l tl+ ...+i n <i TT — = dtx . . . dt n . 

Multiplying and dividing by exp{— a(ti + • • • + £„)} we have 

a n (t) < (Viy 



-at; 



1 e ~as 



(5.13) 



(5.14) 



By choosing a = — we get (5.12). 



T —~=- dt x ... dt n < (Vt) n e a / ^ds <(Vt) n e a (^) n 
[o,i]« H VU lj y/s ] v v /^ / 

(5.15) 
□ 



6. Proof of the Fourier law 

In this section we shall prove Theorem 2.2. The proof relies on Theorem 6.1 below which is 
proved in [3]. Writing A^, n > 1, for the set of all sequences x = (x±, ..,x n ) in with 
distinct entries, we first define the -y-functions 



v x, 



t|// e ) :=E e [ ]{v(xi,t) -p e (xi,t)} 



x E A^, n > 1 



(6.1) 



8=1 



where the process starts with a product measure p e , in particular a single configuration, and 
p e (x,t) is the solution of (2.5) with initial datum p e (x, 0) = //[^(x, 0) = 1]. 

Theorem 6.1. There exist r > and c* > so that the following holds. For any (3* > 

and for any positive integer n there is a constant c n < oo so that for any e > 0, any initial 
product measure p e 

' c n {e-H)- c * n , 



t < e 

sup \v e (x,t\u e )\ < { "\_ « ' „„ 



(6.2) 



a-eA 



Proof of (2.19). Recalling (2.17) and applying Theorem 6.1 we have: 

limjf(t) = -Ul - p(l : t) K ), \im&\t) = -Ul - (1 - p(-ht)) K ) (6.3) 

e— S-U Z e— »0 Z 



and (2.19) follows from (2.14). 
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Proof of (2.18). We can express p(-,t) using the Green function for (2.11) and get 

p(r, t) = P t p(r, 0) + / {P airil (da)p(l, t - s) + P s , r ,_ 1 (^)p(-1, t - s)}, (6.4) 
Jo 

where Ptp(r, 0) = E r 0)l T>t ) with B(t) the standard Brownian motion starting from 

r, and r the hitting time of {—1, 1}; P s ,r,±i(^ s ) are the corresponding hitting time distri- 
butions. Since p'(r,t) := dp(r,t)/dr satisfies the heat equation we can write similarly to 
(6.4) 

p'(r, t) = P t p'(r, 0) + / {¥ s>r>1 (ds)p'{l, t-s) + F s>r> ^{ds)p'(-1, t - s)} (6.5) 
Jo 

with p'(±l,t — s) explicitly given in (2.14). The idea then is to write t) (which is 

defined in (2.16)) in a similar way. We are going to prove that 

= _ * J* | e w (t _ s )P Xie;+ ( rfs ) + eL e) (t - *)P, iej -(<fa)} (6.6) 

with 

Q%\t) = e-^E^N - K,t)] - E e [r,(N - K - 

eL e) (t) = e" 1 (E e [r)(-N + K + 1, t)] - E e [r)(-N + K, t)]) (6.7) 

and 

£ (x + 1, t) = e" 1 ^ (p x+ i{y(t) =y,r>t)- p x { y {t) = yy> t))p*(y, o), (6.8) 

2/ 

where we have used the following notation: y(t) is a random walk with transition kernel 
P^ e \x, y) on Atv, P x its law when y(0) = x, r its first hitting time of {— N + K + 1, N — K} 
and r' the hitting time of {-N + K, N - K - 1}. 

To prove (6.6) we use (2.3), observing that when x ^ J_ U I+, the second term on the r.h.s. 
of (2.3) vanishes, so that 

^E 6 fo(M)] = ^A £ E e [77(x,t)], |x| <N-K. 

This allows to express ^(x, t) := E e [r)(x, t)] = P x [g e (y(fAt),t — f At)] where here P x refers to 
the expectation with respect to a random walk y(-) that starts at x (time running backward 
for y(-)) and f is the first time it reaches {—N + K, N — K}, which corresponds to the (time 
variable) boundary condition. Doing the same for each of the terms in E e [rj(x, t)] — E e [q(x + 
l,t)] we arrive to (6.6). 

We need to compare (6.5) and (6.6) recalling (2.14). By the weak convergence of the random 
walk to the Brownian motion (f)(x + l,t), x = [e _1 r], converges to P t p' (r) and F X!e .±(ds) 
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converges weakly to ¥ rj -±(ds) (x = [e~ 1 r]). Therefore, recalling (2.14), (6.5) and (6.6), the 
proof of (2.18) will follow from: for any t > 



HmeP(t) = j(l -p(l,t) K ), lim6 W (t) = j(l - (1 - P (-l,t)) K ) (6.9) 



which will be proved in the remaining part of this section. As the analysis of 0± (£) are 
similar we shall only prove (6.9) for e { +\t). 

Recalling (2.3) we can write: 

B^(t) = <P e (N - K,t) + ^e,t,y ~ Yl T ^V ( 6 - 10 ) 

yel+ yel- 

where 

<f> e (x,t) e- 1 Y ( p t ie) &y) ~ p t ie) (x -hy))pe(y,o), xeA N ,t>0 (6.11) 

and for y G I±, respectively, 
T eAy := e" 2 J ds (P 5 W(JV - K, y) - P®(N - JT - 1, y)) E e Q(D ± t/(-, t - *))(</)) . (6.12) 

As the analysis of (6.10) will involve several steps, we give first an outline. 

• We shall first prove that <f) e (N — K, t) vanishes as e — > (this will be simple). 

• We will then show that also T e ^,y with y G J_ vanishes as e — )■ 0. This is less simple 
and involves couplings of random walks. 

• The analysis in the previous step is then used to prove that 

lim \T e t y — r* t | =0, for all y G i+, where: (6.13) 

K,t,y ■= e " 2 J ds ( P s e) ( N -K,V)~ P^\N -K-l, y)) 3 -(D +Pe (; t - s))(y). 

• It is then proved that there exist numbers a(h), h = 0, . . . , K — 1, so that 

. K-l 

lim Y K,t,y = \ E a ^ ( X " ^ ( 6 - 14 ) 

ye/+ h=o 

• The final step consists in recognizing that the right hand side of (6.14) is indeed 
equal to j(l — p(l, t) K ). 

By (4.1) 

t) = e" 1 ^ Ql e \x, z) (p e (Mz) 7 0) - Pe (iP N (z - 1), 0)) . (6.15) 
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Recalling that p e (y, 0) = po{ey) and that p' , the derivative of po, is by assumption bounded, 
we then have for any r G [—1,1] 

lim0 e ([e-V],t)= / G t (r,r')(-l) s ^p' ^(r'))dr' =: <f)(r,t), (6.16) 

where S(r') = 1 if r' in [—1,1], ±[3, 5],. . . and = —1 in the complement. By symmetry 
<p(±l,t) = so that 

\im(f) e (x,t) = 0, x = N-K. (6.17) 



By rescaling the time we rewrite T e>t>y as 

»€" 2 t 



|E e [|(D±77(-,t-e 2 s))(y)] (6.1J 



and recall that P^\x,y), which in this proof we denote by ps"'(x,y), is the transition 
probability of a reflected random walk in Ajy with jump intensity 1/2 for each pair of n.n. in 
Ajy. In the sequel we shall also consider the transition probabilities p t (x,y) of the random 
walk on Z + = {0, 1, 2, . . .} with jump intensity 1/2 among nearest neighbors. 



Lemma 6.2. There exists a constant c so that for any h = 0, .., K — 1 and any t 

\p t (K,h)-p t (K+l,h)\ < 



t 3/2- 



The integrals below are well defined: 

POO 

/ {p t (K,h)-p t (K + l,h))}dt=:a(h). 
Jo 



(6.19) 



(6.20) 



Proof. The second statement follows at once from the first, which we now prove with a 
coupling argument. We write 

Pt(K, h) -p t (K + l,h) = £ [l yi (t)=h - l w (t)=ft] , 

where £ is the expectation in a process which couples two simple random walks on Z + , 
denoted by yi(s) and 2/2(5), s G [0,t], with yi(0) = K, 2/2(0) = K + 1. The coupling (whose 
law will be denoted by V) is defined as follows: 2/2(5) moves as the random walk on Z + (i.e. 
with transition probability p s (x,y)) for all s G [0,t]. Let t\ = t/3: in the time interval [0, ti], 
yi(s) copies exactly the jumps of 2/2(5) for all s < min{r, t\}, where r is the first time when 
2/2 jumps to 0. If r < t\ then we set 2/1(5) = 2/2(5) for all s G [r,t]. When r > ti, we let 2/1 
move independently of 2/2 in [ti, t*], where r* is the first time when 2/1 = 2/2 ; and for s > t* 
we set 2/1(5) = 2/2(5). V is evidently a coupling and we have: 



Pt(K, h) - p t (K + l,h) = £ [(l yi (t)=h - l y2 (t)=h)lr*>t] 
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Letting t 2 = 2i/3 and 

g(z u z 2 ) := S [l w (t)=A + l w (t)=h | yi(<2) = z 1 ,y 2 (t 2 ) = z 2 ] 

h(zi,z 2 ) ■= V [t* > t 2 ; | yi(ti) = zi,y 2 (*i) = ^2] , 
the l.h.s. of (6.19) is bounded by: 

and (6.19) follows after recalling that — y 2 {t\) — 1 if r > t±. □ 

Lemma 6.3. There is a constant c so that for any h — 0, .., K — 1, any N and any t 

\pf\N — K,N — ti) — p { t N) (N - K - 1, N - h)\ < T -^. (6.21) 

Moreover for any t < N and c suitable positive constant, 

\pf\xiN - h) -p t (N - x,h)\ <ce~ 5N , x = N-K,N-K-l. (6.22) 

Proof. The same argument used in the proof of Lemma 6.2 proves (6.21). Details are 
omitted. As for (6.22), just notice that pt(N — x, h) is the probability for the random walk 
on {y G Z: y < N} reflected at N and starting at x at time to be at iV — h at time t, 
while Pt(x, N — h) refers to the walk that is also reflected at —N. Letting the two walks 
move together before reaching —N, the difference on the l.h.s is bounded from above by 
the probability of reaching —N by time N, and the estimate follows at once by very simple 
exponential bound on the Poisson clock process (or still using (4.2)). □ 



Lemma 6.4. For any ?/6L 



limr- >y = 0, ye I- (6.23) 

6— >0 ' ' y 



Proof. By (6.18) 

\r e , t , y \<c / \pf\N - K,y) - pi N \N - K - l,y)\ds. (6.24) 
</o 

We bound the probability difference by coupling the two random walks as in the beginning 
of the proof of Lemma 6.2, namely the random walk yi(s) starting at N — K copies the 
jumps of y 2 {s), the one starting at TV — K — 1. Calling the first hitting time of N by y 2 , 
the two random walks become identical after tn- Let r_ be the first hitting time of —N + K 
by y 2 . Thus the contribution to (6.24) comes from the event t_ < t^. Calling V the law of 
the above coupling, S its expectation and T T _ the canonical a-algebra, we have 

|r e ,t,j/| — c ^ [1t_<s1tjv>t_ \£ [lj/i( s )=j/ — l y2 (s)=y\ J~. r _] I] ds. (6.25) 

^0 
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Since tn > r_, yx{rJ) = — N+K+l and ^(t-) = —N+K, the above conditional expectation 
can be bounded using (6.21) (changing x to —x). Thus 



< C 



ds 



(6.26) 



" <s TJV>r -l + (s-r_) 3 / 2 _ 

The r.h.s. of (6.26) involves only the random walk 1/2, and due to the initial conditions we 
are considering (K is fixed), V{j- < tjv) < ce for a positive constant c. Calling m(dt) the 
law of r_ conditioned to r_ < T/v, we may write for y £ I_ 



r e ,t, v | < c'e / ds / m(rftt) 
m(du 



<c'e 



(0,e" 2 t] 



1 + (s - it) 3 / 2 



- <u) 3 / 2 
rfs < c"e. 



proving the lemma 



□ 



Proof of (6.13). We split the integral in (6.18) into s < e~ c * and s > e~ c *, where c* is as in 
Theorem 6.1 (assuming without any loss that c* < 2). For the second one we use (6.21) to 
see that 



< Ce c */ 2 



f ^ ds (pi N \N -K,y)- pi N \N - K - 1, y)) E e ((D ±V (; t - eh))(y)) 
Same estimates hold for r* tJ/ so that using Theorem 6.1 we get 

|pf ) (N -K,y)- p[ N ) (N - K - 1, y) \ds + 2C7e c */ 2 



IF F* I < r( r(2-/?*)c* 

hence (6.13). 

Proof of (6.14). Again by Lemma 6.3 we have 

r 



□ 



(pi N \N -K,y)- pi N \N - K - 1, y)) |(L> +Pe (-,t - e 2 s)){y)ds 



We replace (D + p e (-,t — e s))(y) by (D + p e (-,t))(y), the error being bounded by C\/ -, by 
(5.2). Using again Lemma 6.3 we obtain 

K,t,y-kD + Pe(;t))(y) r (pf\N-K,y)-pi N \N-K-l,y))ds ■ 



By (2.7), 



\im( D + p e (-,t))(y) = (l-p(l,t))p(l,t) N -« ye I, 
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which, by (6.20), proves (6.14). □ 

We are left with the final step, namely to recognize that the right hand side of (6.14) is equal 
to 1 — p(l,t) K . We use conservation of mass, namely from (2.5) it follows that 



N N 



26 £ p e (x,t + r)- PeM) = f i- 1 -J { :\s)+J2 3 2 D ^(y^)) ds 

\x=N-K x=N-K / * \ yel+ J 

(6.27) 

with J+\s) the analogue of 0+\s), namely 

: = e -i ( Pe ( N -K,t)- p e (N - K - 1, t)) . 

Then, analogously to (6.10) 

4\t) = - K,t) + £ r: Ay - J2 K,t, y - 

y ei+ yel- 

The same arguments used for T et v , for y E 7_ shows that lim > T*. = 0, so that using 



(6.14) we get from (6.27) in the limit e — >■ 



2/6/- 



= I E (-^WC 1 - s) ft + (1 - p(l, S ) h ^ cfe (6.28) 



7i=0 

which by the continuity in t gives for any t > 



5 J>(/i)(l - p(l,t))p(l,t) h = £(1 -p(l,f))p(M) fc = 1 -P(M)* (6-29) 



2 

h=0 h=0 



□ 
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